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Monolayer graphene exhibits many spectacu¬ 
lar electronic properties, with superconductiv¬ 
ity being arguably the most notable exception. 
It was theoretically proposed that supercon¬ 
ductivity might be induced by enhancing the 
electron-phonon coupling through the decora¬ 
tion of graphene with an alkali adatom super¬ 
lattice [Profeta et aL Nat, Phys, 8, 131-134 
(2012)]. While experiments have indeed demon¬ 
strated an adatom-induced enhancement of the 
electron-phonon coupling, superconductivity has 
never been observed. Using angle-resolved pho¬ 
toemission spectroscopy (ARPES) we show that 
lithium deposited on graphene at low temper¬ 
ature strongly modifies the phonon density of 
states, leading to an enhancement of the electron- 
phonon coupling of up to A 0.58. On part 
of the graphene-derived 7r*-band Fermi surface, 
we then observe the opening of a A 0.9 meV 
temperature-dependent pairing gap. This result 
suggests for the first time, to our knowledge, that 
Li-decorated monolayer graphene is indeed super¬ 
conducting with Tc^5.9K. 

While not observed in pure bulk graphite, supercon¬ 
ductivity occurs in certain graphite intercalated com¬ 
pounds (GICs), with Tc of up to 11.5K in the case of 
CaCe m E]- The origin of superconductivity in these 
materials has been identified in the enhancement of 
electron-phonon coupling induced by the intercalant lay¬ 
ers lami. The observation of a superconducting gap on 
the graphitic tt* bands in bulk CaCe [5] suggests that re¬ 
alizing superconductivity in monolayer graphene might 
be a real possibility. This has indeed attracted intense 
theoretical and experimental efforts ISHI2- In partic¬ 
ular, recent density-functional theory calculations have 
suggested that, analogous to the case of intercalated bulk 
graphite, superconductivity can be induced in monolayer 
graphene through the adsorption of certain alkali metals 
0 . 

Although the Li-based GIG - bulk LiGe - is not known 
to be superconducting, Li decorated graphene emerges as 
a particularly interesting case with a predicted supercon¬ 
ducting Tc of up to 8.IK [8]. The proposed mechanism 
for this enhancement of Tc is the removal of the confin¬ 
ing potential of the graphite Ge layers, which changes 


both the occupancy of the Li 2s band (or the ionization 
of the Li) and its position with respect to the graphene 
layer. This ultimately leads to an increase of the electron- 
phonon coupling constant from A = 0.33 to 0.61, in going 
from bulk to monolayer LiGe- It has been argued that 
the LiGe monolayer should exhibit the largest values of 
both A and Tc among all alkali-metal-G6 superlattices 
[8]. Nevertheless, while there is thorough experimental 
evidence for adatom-enhanced electron-phonon coupling 
in graphene 13 HD ig, superconductivity has not yet 
been observed in decorated monolayer graphene. 

ARPES measurements of the electronic dispersion of 
pristine and Li-decorated graphene at 8K, character¬ 
ized by the distinctive Dirac cones at the corners of 
the hexagonal Brillouin zone [Fig.j^e)], are shown in 
Fig. 0 a) and (b). Li adatoms electron-dope the graphene 
sheet via charge-transfer doping, leading to a shift of the 
Dirac point to higher binding energies. As evidenced 
by the evolution of the graphene sheet carrier density 
in Fig.j^f), this trend begins to saturate after several 
minutes of Li deposition. Goncomitantly, we observe 
the emergence of new spectral weight at the Brillouin 
zone centre [see Fig.j^e) and the comparison of the T- 
point ARPES dispersion for pristine and 10 minute Li- 
decorated graphene in Fig.[^c,d)]. The origin of this 
spectral weight is probably the Li-25 band expected for 
this system [8], superimposed with the folded graphene 
bands due to a Li superstructure, as observed in Li and 
Ga bulk GIG systems [5l[T4]. This spectral weight, which 
disappears above ^ 50 K and is not recovered on subse¬ 
quent cooling (see SI Appendix), is associated with the 
strong enhancement of electron-phonon coupling to be 
discussed later. 

Next we use high-resolution, low-temperature ARPES 
to search for the opening of a temperature-dependent 
pairing gap along the 7r*-band Fermi surface, as a direct 
spectroscopic signature of the realization of a supercon¬ 
ducting state in monolayer LiGe- To increase our experi¬ 
mental sensitivity, as illustrated in Fig. ^a) and following 
the approach introduced for FeAs [17 and cuprate [18] 
superconductors, we perform an analysis of ARPES en¬ 
ergy distribution curves (EDG) integrated in dk along a 
one-dimensional momentum-space cut perpendicular to 
the Fermi surface. This also provides the added benefit 
that the integrated EDGs can be modelled in terms of 
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FIG. 1. Charge-transfer doping of graphene by lithium adatoms, (a) Dirac-cone dispersion measured by ARPES at 8K 
on pristine graphene and (b) after 3 minutes of Li evaporation, along the iF-point momentum cut indicated by the white line 
in the Fermi surface plot in (e). The Dirac cone Fermi surface was measured at this specific K point, and then replicated at the 
other K points by symmetry (note that high-symmetry points are here defined for the Brillouin zone of pristine graphene and 
not of a/3 X a/3R30° reconstructed Li-graphene, which is instead the notation in Ref. [8}. The point at which the spectroscopic 
gap is studied is indicated by the shaded white circle. The Dirac point, already located below Ep on pristine graphene due 
to the charge-transfer from the SiC substrate (a), further shifts to higher energies with Li evaporation (b). The presence of 
a single well-defined Dirac cone indicates a macroscopically uniform Li-induced doping. While no bands are present at the 
T-point on pristine graphene (c), spectral weight is detected on 10-minute Li-decorated graphene in (d) and (e). As illustrated 
in the 8K sheet carrier density plot versus Li deposition time in (f), which accounts for the filling of the tt* Fermi surface, the 
spectral weight at T is observed for charge densities n2D^9x lO^^cm”^ (but completely disappears if the sample temperature 
is raised above ~50K, and is not recovered on subsequent cooling; see also SI Appendix). 


a simple Dynes gap function [19] multiplied by a linear 
background and the Fermi-Dirac distribution, all con¬ 
volved with a Gaussian resolution function (see Methods 
and in particular Eq.|^. As shown in Fig.j^a) and espe¬ 
cially l^b) for data from the k-space location indicated 
by the white circle in Figs. [^e) and [^e), a temperature 
dependence characteristic of the opening of a pairing gap 
can be observed near Ep- The leading edge midpoint of 
the Li-graphene spectra moves away from Ep [Fig.j^b)] 
in cooling from 15 to 3.5 K, at variance with the case 
of An spectra crossing precisely dX Ep according to the 
Fermi-Dirac distribution [Fig.j^d)]. When fit to Eq.|^ 
this returns a gap value A = 0.9±0.2meV at 3.5 K (with 
F 0.09 meV [20]). Given its small value compared to 
the experimental resolution, the gap opening is best vi¬ 
sualized in the symmetrized data in Fig.j^c), which min¬ 
imizes the effects of the Fermi function. Finally, we note 
that the gap appears to be anisotropic, and is either ab¬ 
sent or below our detection limit along the K—M direc¬ 
tion (see Fig. 4 in SI Appendix). 

The detection of a temperature-dependent anisotropic 
gap at the Fermi level with a leading-edge profile de¬ 
scribed by the Dynes function - with its asymmetry 
about Ep and associated transfer of spectral weight to 
just below the gap edge - is suggestive of a supercon¬ 


ducting pairing gap. The phenomenology would in fact 
be very different in the case of a Coulomb gap, typi¬ 
cally observed in disordered semiconductors [24U26] , due 
to the combination of disorder with long-range Coulomb 
interactions. This would lead to a rigid shift of the spec¬ 
tra leading edge, isotropic in momentum, and result in a 
vanishing of the momentum-integrated density of states 
at E^. 

Similarly, the observed gap is unlikely to have a charge 
density wave origin, since the gap is tied to the Eermi en¬ 
ergy as opposed to a particular high-symmetry wavevec- 
tor (the latter might occur at the M points, when 
graphene is doped all the way to the Van Hove singulari¬ 
ties resulting in a highly-nested hexagonal Eermi surface, 
or at the K points in the case of a a/3 x a/3 R 30° recon¬ 
struction leading to a Dirac-point gap). Einally, we note 
that these measurements do not allow us to speculate on 
the precise symmetry of the gap along a single Dirac-cone 
Eermi surface, nor on the relative phases of the gap on 
the six disconnected Eermi pockets. As such, our results 
do not rule out any of the recent proposals for a possible 
unconventional superconducting order parameter (see for 
example Refs.[9l [271 and [28]) . 

To further explore the nature of the gap observed on Li- 
decorated graphene (and also demonstrate our ability to 
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FIG. 2. Spectroscopic observation of a pairing gap in Li-decorated graphene, (a) Dirac dispersion from 10-minute 
Li-decorated graphene measured at 15 and 3.5 K, at the k-space location indicated by the white circle in Figs. [^e) andlUe); 
the temperature dependence is here evaluated for EDCs integrated in the 0.1 momentum region about kp shown by the 
white box (bottom), with the only changes occurring near Ep (top). While An spectra (d) cross at Ef as described by the 
Fermi-Dirac distribution, the crossing point of the Li-graphene spectra (b) is shifted away from Ep (cyan dashed line), due to 
the pull-back of the leading edge at 3.5 K. A fit to the Dynes gap equation (see Methods) yields a gap A~0.9meV at 3.5 K (and 
OmeV at 15 K). The superconducting gap opening is best visualized in the symmetrized data in (c), i.e. by taking 
which minimizes the effects of the Fermi function even in the case of finite energy and momentum resolutions [151116j [blue and 
red symbols in (c) represent the smoothed data, while the light shading gives the root-mean-square deviation of the raw data]. 
The qualitatively similar behaviour observed on polycrystalline niobium - and returning a superconducting gap A~1.4meV - 
is shown in (e) and (f). 


resolve a gap of the order of 1 meV), in Fig.[^e,f) we show 
as a bench-mark comparison the analogous results from 
a bulk, polycrystalline niobium sample - a known BCS 
superconductor with Tc — 9.2K. The Dynes fit of the 
integrated EDCs Fermi edge in Fig.j^e) determines the 
gap to be A = 1.4±0.2 meV (with F:^0.14meV [20]), in 
excellent agreement with reported values [29]. Although 
the leading edge shift [Fig.j^e)] and the dip in the sym¬ 
metrized spectra [Fig.j^f)] are more pronounced than for 
Li-graphene owing to the larger gap, the behaviour is 
qualitatively very similar. This provides additional sup¬ 
port to the superconducting origin of the temperature- 
dependent gap observed in Li-decorated graphene. 

If this is indeed a superconducting gap, the responsi¬ 
ble mechanism may likely be electron-phonon coupling, 
as predicted by theory for monolayer Li-graphene [8] and 


also seen experimentally for bulk GIC CaCe [5] . In direct 
support of this scenario, we present a detailed analysis of 
the graphene tt* bands in Fig.|^ demonstrating that the 
Li-induced enhancement of the electron-phonon coupling 
is indeed sufficient to stabilize a low-temperature super¬ 
conducting state. Graphene doped with alkali adatoms 
always shows a strong kink in the tt* band dispersion 
at a binding energy of about 160 meV m- For the Li- 
graphene studied here, this is seen in the momentum- 
distribution curve (MDC) dispersions and correspond¬ 
ing real part of the self-energy E' in Fig.[^b-d). This 
structure stems from the coupling to carbon in-plane 
{Cxy) phonons [4] [8]. Despite the apparent strength of 
this kink, the interaction with these phonon modes con¬ 
tributes little to the overall coupling parameter due to 
their high energy (note that uj appears as a weighting 
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FIG. 3. Analysis of electron-phonon coupling in Li-decorated graphene, (a) Dirac dispersion from 3-minute 
Li-decorated graphene, along the /c-space cut indicated in the Fermi surface plot in (e), exhibiting kink anomalies due to 
electron-phonon coupling (white line: MDC dispersion), (b-d) MDC dispersion and bare-band obtained from the self-consistent 
Kramers-Kronig bare-band fitting (KKBF) routine [21] [22], for several Li coverages (see Methods and SI Appendix); the real 
part of the self-energy E' is shown in the side panels (orange: E' from KKBF routine analysis; black: E' corresponding to the 
Eliashberg function presented below), (f-h) Eliashberg function a^F((jj) from the integral inversion of E'(ct;) [23], and electron- 
phonon coupling constant \ — 2 Jduj F{uj)/u (see Methods and SI Appendix); in (h) the theoretical result from Ref. [8] for 
a LiCe monolayer are also shown (gray shading), (i) Experimentally-determined contribution to the total electron-phonon 
coupling (black symbols) from: phonon modes in the energy range 100-250meV (blue shading, white symbols), and 0-100meV 
(orange shading); the coupling of low-energy modes strongly increases with Li coverage. 


factor in the integral calculation of A - see Methods). As 
illustrated by the white symbols in Fig.j^i), the contri¬ 
bution to A from these high-energy (100-200 meV) modes 
is determined to be 0.14±0.05, and it remains approxi¬ 
mately constant for all Li coverages studied here. This 
value is, however, too small to stabilize a superconduct¬ 
ing state in this system min]. 

With increasing Li coverage and the appearance of the 
spectral weight at F, significant modifications to the low- 
energy part of the dispersion (< 100 meV) become ap¬ 
parent [Fig.[^b-d)]. With 10 minutes of Li deposition 
[Fig-id)], an additional kink is visible at a binding en¬ 
ergy of approximately 30 meV, along with the associated 
peak in the real part of the self-energy E'. The extracted 
(see Methods) Eliashberg functions and energy-resolved 
A(ci;) in Fig.f-h) show that, at high Li coverage, phonon 
modes at energies below 60 meV are coupling strongly to 
the graphene electronic excitations. The phonon modes 
in this energy range are of Li in-plane (Li^,^) and C out- 
of-plane (Cz) character [HIS]. This is in agreement with 
predictions [S], as shown by the direct comparison be¬ 
tween theory and experiment in Fig.j^h) [30]. As for the 
total electron-phonon coupling A for each coverage [black 
symbols in Fig.j^i)], our values measured on the 7r*-band 
Fermi surface at an intermediate location between T—K 
and K—M directions [Fig.|^e)] provide an effective esti¬ 
mate for the momentum-averaged coupling strength m- 


Remarkably, the value A = 0.58 ± 0.05 observed at the 
highest Li coverage [Fig.j^i)] is comparable with A = 0.61 
predicted for monolayer LiCe [8] as well as A 0.58 ob¬ 
served for bulk CaCe [32] - it is thus large enough for 
inducing superconductivity in Li-decorated graphene. It 
is also significantly larger than the momentum-averaged 
results previously reported for both Li and Ca deposition 
on monolayer graphene (A 0.22 and 0.28, respectively 
my We note that achieving such a large A value is crit¬ 
ically dependent on the presence of the spectral weight 
observed at F when Li is deposited on graphene at low 
temperatures, presumably forming an ordered structure 
on the surface and not intercalating. As shown in the SI 
Appendix, we find A = 0.13 ± 0.05 after the same sample 
is annealed at 60 K for several minutes, destroying the Li 
order and associated F spectral weight. 

Taken together, our ARPES study of Li-decorated 
monolayer graphene provides the first evidence for the 
presence of a temperature-dependent pairing gap on part 
of the graphene-derived tt* Fermi surface. The detailed 
evolution of the density of states at the gap edge, as well 
as the phenomenology analogous to the one of known su¬ 
perconductors such as Nb - as well as CaCe and NbSe 2 , 
which also show a similarly anisotropic gap around the K 
point [33H37] - indicate that the pairing gap observed at 
3.5 K in graphene is most likely associated with supercon¬ 
ductivity. Based on the BCS gap equation, A = 3.5 Tc, 
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this suggests that Li-decorated graphene is superconduct¬ 
ing with Tc ^ 5.9 K, remarkably close to the value of 8.1 K 
found in density-functional theory calculations [8]. This 
constitutes the first experimental realization of supercon¬ 
ductivity in graphene - the most prominent electronic 
phenomenon still missing among the remarkable proper¬ 
ties of this single layer of carbon atoms. 
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METHODS 

Sample preparation. Epitaxial graphene monolayers 
with a carbon buffer layer were grown under argon at¬ 
mosphere on hydrogen-etched 6H-SiC(0001) substrates, 
as described in Ref.|3Hl The samples were annealed at 
500° C and 8 X10 Torr for 1 hour, immediately prior 

to the ARPES measurements. Lithium adatoms were 
deposited from a commercial SAES alkali metal source, 
with the graphene samples held at a temperature of 8K. 
Bulk Nb polycrystalline samples, with Tc = 9.2 K, were 
cleaved in the ARPES chamber prior to the experiments. 
ARPES experiments. The measurements were per¬ 
formed at UBC with 5-polarized 21.2 eV photons on an 
ARPES spectrometer equipped with a SPECS Phoibos 
150 hemispherical analyzer, a SPECS UVS300 monochro- 
matized gas discharge lamp, and a 6-axes cryogenic ma¬ 
nipulator that allows controlling the sample temperature 
between 300 and 3.5K, with accuracy TO. IK. Band and 
Fermi surface mapping, as well as the study of electron- 
phonon coupling, were performed at 8 K with energy and 
angular resolution set to 15 meV and 0.01 A“^, respec¬ 
tively. For the measurements of the superconducting 
gaps, energy and angular resolution were set to 6 meV 
and 0.01 A“^, while the sample temperature was varied 
between 3.5 and 15 K. During the ARPES measurement 
the chamber pressure was better than 4x10“^^ Torr. 
Electron-phonon coupling analysis. The spectral 
function A(k, uj) measured by ARPES [39] provides infor¬ 
mation on both the single-particle electronic dispersion 
e\^ (the so-called ‘bare-band’) as well as the quasiparticle 


self-energy E(k, cj) = E'(k, co’)+'iE"(k, cj), whose real and 
imaginary parts account for the renormalization of elec¬ 
tron energy and lifetime due to many-body interactions, 
including electron-phonon coupling. By fitting with a 
Lorentzian and a constant background the ARPES inten¬ 
sity profiles at constant energy u; = a), known as momen¬ 
tum distribution curves (MDGs), one obtains the MDC 
dispersion defined by the peak maximum km [plotted in 
Fig-i a-d)], as well as the corresponding half-width half¬ 
maximum (HWHM) Akm- The real and imaginary parts 
of the self-energy can then be defined as: 

Y.l = -^kraVl^, ( 1 ) 

(where is the bare-band velocity). To extract the 
self-energy and dispersion without any a priori knowledge 
of the bare-band, we use the self-consistent Kramers- 
Kronig bare-band fitting (KKBF) routine from Ref. [2] 
and[2^ (see also SI Appendix). 

As for the dimensionless /c-resolved electron-phonon 
coupling constant discussed in the paper and in particu¬ 
lar in Fig. if- i), this is formally defined as [40] : 


nUJ 

Ak(ci;) = 2 / duo' 

Jo 


a^F{k,u') 


( 2 ) 


where (a^F(k, cj) is the Eliashberg function, i.e. the 
phonon density of states weighted by the electron-phonon 
coupling strength [40]. The latter is related to the real 
part of the self-energy E'(k, cj) via the integral relation: 





duj' a^F{k,uj')K 


/ UJ cj' \ 


( 3 ) 


where K{y,y') = j^^dxf{x - y)2y'l{x‘^ - y'‘^) and 
f{x — y) is the Fermi-Dirac distribution. The momen¬ 
tum resolved a^Fik^uj) function plotted in Fig.|^f-h) 
can then be extracted from the real part of the self¬ 
energy E'(k, cj) probed by ARPES via the integral in¬ 
version procedure described in Ref.|23l Ultimately, by 
means of Eq.|^ this also allows the calculation of the 
electron-phonon coupling constant shown in Fig.[^ (see 
also SI Appendix). 

Superconducting gap fitting. As discussed in detail 
for the case of FeAs and cuprate superconductors in 
Ref.[T7]and |T8l respectively, when evaluating the opening 
of a superconducting gap based on the EDCs integrated 
in dk along a one-dimensional momentum-space cut per¬ 
pendicular to the Fermi surface (such as those presented 
in Fig.[^, one can make use of the following formula: 




f{uj,T) (aFbuj) 


Re 


UJ — iV 


v/(cj-ir)2-A2 


^Ru. 


( 4 ) 

This corresponds to the Dynes gap function (i.e., the 
Bardeen-Cooper-Schrieffer density of states with a su¬ 
perconducting gap A, broadened by the pair-breaking 
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scattering rate rm), multiplied by a linear background 
(with parameters a and b) and by the Fermi-Dirac distri¬ 
bution all convolved with a Gaussian function 

accounting for the experimental energy resolution 
(owing to the integration of the ARPES intensity in dk, 
this analysis is unaffected by momentum resolution in])- 
The high and low temperature ARPES data in Eig.j^are 
fitted simultaneously using the above equation, accord¬ 
ing to the following additional considerations: since the 
integrated EDCs were observed not to change outside of 
the gap region, the linear background is constrained to 
be the same for the above and below Tc measurements; 
the temperatures are fixed to the known measured values 
(with accuracy iO.l K); the Eermi energy Ep and energy 
resolution are determined from fitting the high temper¬ 
ature data from either Nb and Li-graphene during each 
measurement, and independently verified from measure¬ 
ments on polycrystalline gold. 
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